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ABSTRACT 

We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of 
blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our 
measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble 
Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope ob- 
servations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, 
suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In 
the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, 
making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galac- 
tocentric distances. If we assume a power-law profile r~ a for the 2-D projected surface density BHB 
distribution, we obtain a high-quality fit with a 2-D power-law index of a = 2.6^0 2 outside of 3 
kpc, which flattens to a <1.2 inside of 3 kpc. This slope is consistent with previous measurements 
but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal 
symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass 
of 2.llg'4 x 10 9 Mq. These properties are comparable with both simulations of stellar halo formation 
formed by satellite disruption alone, and with simulations that include some in situ formation of halo 
stars. 

Subject headings: galaxies: individual (M31) — galaxies: stellar populations — galaxies: evolution 



1. INTRODUCTION 



outsk irt s of many nearby galaxie s ( e.g.. iMalin fe Hadlev 



The diffuse envelope of stars that appears to sur- 
round most or all Milky Way-mass galaxies (massive 
galaxies, hereafter) encodes a wealth of information 
about how such ga l axies assemble d and d eveloped (e.g. 
Eggen et all \l96% iSearle fe Zinnl fl97l IBullock et all 



20011: IBullock fc Johnst on 2001 [Zolotov et all 120091: 



Cooper et alll2010l ). A major fraction of the stellar mass 



in such envelopes is expected to be debris from tidally 
disrupted dwarf galaxies embedded in the dark matter 
halos. These disrupted galaxies are progressively incor- 
porated into the larger potential well of the main galaxy 
during its hierarchical, bottom-up assembly Accord- 
ingly, widespread effort has been put into studying mas- 
sive galaxy halos both observationally and theoretically. 

Over the past few decades, major observational 
progress has been made in our understanding of the 
build-up of halos through discovery of many dwarf 
galax ie s around the Milky Way (e.g.. Ilbata et all 
[1991: IWillman et all 120051 IBelokurov et all 120071) . 
low surface b rightness structure s in the Milk y Way 
stellar halo (Ilbata et~aTI 119951: lYannv et all 120001 : 
iMaiewski et all 120031: iBell et all l2008l ).and around the 
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IMartinez-Delgado et~aTI l200cl iBailin et al 
iRadburn-Smith et all 120 111), including M31 (e.g. 



Ibata et al.ll2001blGilbert et al. 2007: McConnachie et all 
20091: IGilbert et alJl2009aTT The amount and character 
of such substructure at r > 10 kpc is both qualitatively 
and quantitatively consistent with cosmologically- 
motivated simulations in which stellar halos arc 
built u p through the acc r etion o f dwarf g a laxies 
alone (iBullock fc Johnston] 12001 IBell et all l200l 
IMartinez-Delgado et al .1120101: IXueet al.H2011[) . 

At radii less than ~20 kpc, it is expected that some 
fraction of stellar halo mass may originate in the main 
potential well (often called in situ), either formed early 
in the formation of the galaxy or kicked up from 
the disk at later times through tidal interactions (e.g. 
Kazantzidis et all 120081: iZolotov et~aTl 120091 : iFont et all 



20111 : iBrook et all 120121 ). Studies attempting to disen 
tangle the fraction formed in situ 1-2 are challenging both 
observationally and theoretically for a variety of reasons. 

Observationally, the stellar halo has both low surface 
brightness and relatively low stellar mass, and so is eas- 
ily overwhelmed by higher surface brightness components 
of a galaxy, making accurate measurements challenging. 
A further (related) difficulty is that the separation of a 
halo from a bulge component (where the bulge in great 
part is expected to have formed the bulk of its stars 
in the main galaxy potential well) has been extremely 
diffic ult because disk and bulge dominate the stellar 
halo (jDorman et all l2012t iCourteau et~aT1 12011L here- 
after Cll), leading a nu mber of groups to consider both 
components together (e.g.llrwin et aTi ; 2005; K alirai et all 



1200a iSaraiedini et alll2012l: IFont et al. 
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A related theoretical problem is that in numerical sim- 
ulations relatively rare modes of star formation at rel- 
atively high temperatures and low densities, which are 
not yet well-characterized, may contribute disproportion- 
ately to low-mass diffuse stellar halos (see discussions in 
IZolotov et aIJ200ftlMunshi et al.l2012t) . The correct star 
formation prescription for the halo gas is still uncertain, 
making the simulated stellar halo properties less reliable. 
In what follows, we attempt to give a fresh perspective to 
this issue by focusing on the possibility of using blue hor- 
izontal branch (BHB) stars in M31 as a probe of M31's 
inner stellar halo. 

1.1. Current M31 Stellar Halo Measurements 

Fits to the surface brightness profile of M31 do not con- 
strain the shape of the inner halo. The surface brightness 
profile is dominated by the red giant branch (RGB) stars 
in M31. These studies fit all three galaxy components 
(disk, halo, and bulge) simultaneously, but must assume 
a halo model inside of — 10 kpc where all 3 components 
contribute significantly. Not only are the fits to the inner 
halo profile model dependent, the fits by Cll show clear 
degeneracies between the different components, making 
extrapolating the halo profile inwards from 10 kpc unre- 
liable. Thus, the modeling of the inner halo can result in 
significant variations of the stellar halo mass fraction at 
these inner radii. 

Kinematic decompositions of the halo outside of 
10 kpc, performed with spectra of individual RGB stars, 
currently suggest a 2-D projected surface density power- 
law index of ——2 (K. Gi lbert et al. ApJ, submitted; 
iGuhathakurta et "all 120051 hereafter G05). However, 
M31's bulge and halo are both dynamicall y hot and ex- 
hibit metallicity and/or a ge gradients (e.g. iKalirai et al.l 
120061: IGilbert et al.l [20071 ) . making it extremely difficult 
to disentangle M31's bulge from M31's halo closer to the 
galaxy center, even with kinematics of individual red gi- 
ant branch stars. 

Furthermore, a case can be made that it is danger- 
ous to extrapolate the outer halo stellar profiles in- 
wards to small galactocentric distances in M31, or in 
any other galaxy. Inside of —35 kpc, the stellar popula- 
tions of the M31 halo are substantially more metal-rich 
than those of the Milky Way and have a range of ages 
(IDurrell et al.lH99l [200ll l2Cjoi IBrow n et al.ll2006l [2007L 
l2008HRichardson et al.ll2008j ). In addition, there are de- 
viations from a single power law in the surf ace brightness 
profil e of the outer stellar envelope of M31 (jKalirai et al.l 

[Mil- 
Taken together, these differences have led some work- 
ers to suggest that at radii < 35 kpc one should think 
of the extended stellar envelope as an extended bulge 
instead of halo, or one should simply lump both com- 
ponents together into a spheroid. In addition, simula- 
tions of stellar halo formation, both those that are com- 
posed only of tidal debris from dwarf disruptions and 
those that include an in situ component, show changes 
in power law profiles with radius, freque ntly flattening 
in profile towards the smallest radii (e.g. IZolotov et al.l 
[2001 ICooper et al.|[2010 l : iFont et al.ll201lD . Without re- 
liable constraints on the shape of the inner halo, com- 
parisons of the observed halo to simulations or estimates 
the stellar halo mass are hampered. 
In short, there arc currently no straight-forward tools 



for mapping the halo profile in the inner galaxy. Thus 
any current knowledge about stellar halo profiles is lim- 
ited to large radii. 

1.2. Using BHB Stars to Trace the Halo 

Our approach to recovering the inner structure of the 
halo is to identify a tracer unique to the stellar popula- 
tion of the halo. More specifically, we attempt to define 
the stellar halo as the population that is responsible for 
the presence of field blue horizontal branch (BHB) stars. 
We then check the consistency between our results and 
previous halo studies. We attempt this new approach in 
order to avoid the difficulties of more direct kinematic or 
photometry measurements. 

The halo measurements of previous studies have been 
limited to the outskirts of the gala xy to ensure they 
are measuring a halo population (e.g . IBrown et al.1120081 : 
IRichardson et al.l 120081: IKalirai et al.H2006t) . In these re- 
gions, previous works have found a total stellar popu- 
lation that is poorly fit with a single power-law profile, 
motivating a description of the outer parts of M31 in 
terms of a stellar halo (often taken as synonymous with 
a power law component), an extended bulge (a Sersic in- 
dex component pr ominent out to large radius), and an 
extended disk (e.g. Ilbata et al]|2005b IKalirai et al"1l2006l 
Cll). 

In contrast to current surface brightness and spec- 
troscopic techniques, the BHB may provide a way to 
trace the halo profile to small radii, breaking degenera- 
cies between the disk, bulge, and halo galaxy compo- 
nents. The BHB is the part of the horizontal branch just 
blueward of the RR Lyrae instability strip. It is made 
of low-mass stars in the phase of central helium burn- 
ing, with effective temperatures from —7200-— 40000 K. 
The BHB is routinely subdivided by temperature into 
three parts, the HBA (-7200 K<T eff <12000 K), HBB 
(12000 K<T off <20000 K), and E HB (T7ff>2000 K) (for 
a full review on HB evolution see [Catalan 200j| and ref- 
erences therein). In this work we include only the HBA, 
which is the section typically associated with metal-poor 
populations, when referring to the BHB. For a typical 
old and metal-poor halo environment, these BHB stars 
will spend their —100 Myr lifetime confined within a hor- 
izontal strip just 0.2 mag wide in Mboi, so that the BHB 
location can be approximated by a simple zero-age hor- 
izontal branch (ZAHB) sequence, as the one illustrated 
in Figure [T] 

Although BHB stars are typically associated with low 
metallicity populations, it is also possible to make hot HB 
stars at high mctallicitics. Enhanced mass loss, and pos- 
sibly high helium content, generate a second formation 
channel; however this metal-rich channel is more likely to 
produce extremely hot HB stars, not normal HBA stars, 
as w itnessed by the me tal-rich open cluster NGC6791 
(e.g. IKalirai et al.l l2007j) . and as demonstrated by the 
very low formatio n efficiency of RR Lyrae among metal- 
rich populations (|Lavdenlll995|) . BHB stars instead, are 
likely to have metallicitie s similar to the RR Lyrae , i.e. 

with a mean [Fc/H] 1 (jPietrukowicz et al.l [20121) and 

just a very minor tail of stars extending to high metal- 
licities (see also Figure [SJ . For this reason, these stars 
have been successfu lly used as probes of the Galactic halo 
stellar population (|Preston et al.l fT991h . As a result, a 
prominent BHB component is unlikely to be formed by 
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the higher metallicity disk and bulge populations. We 
note the possibility of a contribution from the thick disk, 
which has recently been measured to have a metallicity 
similar to that of the halo (jCollins et a l. 2011), but argue 
against this possibility on the basis of previous M31 disk 
pop ulation studies and k inematic arguments (see § 13.7.21 
also lKinman et aT1l2009l who find a very low fraction of 
local BHB stars originating in the thick disk). 

There i s clearly a BHB com ponent in the M31 stel- 
lar halo (jvan den Berghl [T991D . a s easily seen in deep 
HST CMDs in the literature (e.g. IBellazzini et al.l 120031: 
IBrown et al.l [200 31 . Furthermore, RR Lyrae stars have 
been cataloged in M 31 far from the central bulge (e.g., 
iBernard etaLl 12012ft . We note that the BHB stars 
we are using to probe the halo arc distinct from ex- 
treme horizontal branch (EHB) stars, which are much 
hotter that the traditional BHB. These EHB stars are 
com mon in metal-rich pop ulations including the M31 
(e.g. iRosenfield et al.l I2012D. and Galactic bulges (e.g. 
lO'ConnelllllOQa iBusso et al.ll2005D and therefore would 
not be appropriate tracers of the halo. Furthermore, 
EHB stars are quite faint in the optical, as their spectral 
energy distribution peaks in the far-UV, making them 
fall below the completeness limit in our ACS data. 

Herein we use BHB stars to trace the M31 halo compo- 
nent to small galactocentric distances, where the overall 
fraction of halo stars is low and difficult to measure any 
other way. We isolate the BHB in the M31 field popu- 
lation and measure its surface density to galactocentric 
distances as small as 1.6 kpc from the center of M31. 
Our measurement of the halo profile provides the means 
to break degeneracies with other components, as well as 
to estimate the total halo stellar mass. In § 2 we de- 
scribe our data set and measurement techniques. In § 3 
we present the resulting BHB density profile and argue 
that it is reasonable to assume that they are all halo 
members. We then quantify our constraints on the halo 
profile and compare these to currently-available decom- 
position measurements and simulations. Finally, in § 5 
we provide a brief summary of our results and their place 
in our understanding of the structu re of M31. 

We assume a distance of 780 kpc (jStanek fc Garnavichl 
[1991 for all conversions from angular distances to kilo- 
parsecs. All power-law indexes refer to 2-D projected 
surface density profiles. We make no corrections for in- 
clination, as there is no clear evidence in the literature 
that the M31 stellar halo is non-spherical, we have no 
way of knowing if an inclination correction is needed. If 
the halo is indeed significantly flattened and inclined, it 
would affect the magnitude ranges adopted for our BHB 
selection in the outer halo fields by ^0.1 mag. Such an 
offset would not significantly impact our surface density 
measurements; however, flattening would affect our mass 
estimate, as discussed in § 13.41 

2. DATA ACQUISITION AND ANALYSIS 

2.1. Panchromatic Hubble Andromeda Treasury Data 

As part of the PHAT survey (|Dalcanton et al.ll2012l ). 
we obtained HST/ ACS data in F475W and F814W cov- 
ering a large fraction of the M31 disk. All of the PHAT 
data for this paper were acqui red and analyzed as par t 
of the survey, as detailed in iDalcanton et al.l (|2012l ). 
In short, photometry was performed using the package 



DOLPHOT ([Dolphin! 1200(1 . The output was filtered to 
reject non-poin t sources and low-quali ty measurements, 
as detailed in IDalcanton et al.l (|2012f) . Artificial star 
tests were performed on each field to allow measurement 
of completeness as a function of color and magnitude. 

In studying the survey photometry, we noticed a blue 
feature extending across the vertical main-sequence in a 
few of our survey fields (see Figure [1]). Further exami- 
nation of the areas exhibiting this feature revealed that 
they all fell within a small corner of the survey, shown in 
Figure [2j This corner lies at the farthest point that the 
survey reaches along the M31 minor axis, and thus has 
the highest ratio of halo to disk stars. 

2.1.1. Counting BHB Stars 

We isolated the BHB component by taking a vertical 
cut in the CMD at 0.1<F475W-F814W<0.5 (see dashed 
vertical lines in Figure [lj, and plotting F475W and 
F814W magnitude histograms. Assuming £'s_y=0.1, 
this intrinsic color range (0.1<F475W-F814W<0.5 after 
correcting for reddening) corresponds to 7200K-9600K. 
Our major axis disk field luminosity functions using this 
color selection was well-fit with a simple straight line, 
suggesting that the underlying main sequence population 
has a simple linear luminosity function. Wc therefore fit 
these luminosity functions, which clearly contained a sig- 
nificant enhancement at the magnitude expected for the 
BHB, with a straight line plus a Gaussian; the latter 
component was then used as a measure of the magnitude 
and strength of the BHB feature at this color. We chose 
our color range based on the data in Figure [T] At bluer 
colors, the BHB sequence extends vertically in our filter 
set, blending completely with the main sequence. Red- 
ward of our BHB selection window, the RR Lyrae in- 
stability strip spreads the feature, also making the BHB 
difficult to isolate. In F814W, we found the peak magni- 
tude to be consistent across all fields where we measured 
it; however, in F475W, the peak magnitude was slightly 
less consistent, owing to small amounts of dust redden- 
ing. Therefore we show the less dust-affected F814W fits 
in Figure [3] Finally, we attempted rotated LFs that were 
orthogonal to the BHB in this color range. We found in 
these cases the background LF became more complicated 
as did comparisons across different filter sets. In the end, 
the simple F814W luminosity function provided the most 
reliable measurements across our full sample. 

The line component of our fits represents the up- 
per main-sequence stars that also occupy this color- 
magnitude region. The number of BHB stars was ob- 
tained by subtracting the line component away from the 
total, thus removing the upper main-sequence contami- 
nation from our BHB sample. 

We performed checks to ensure that our sample was 
not strongly affected by the presence of dust in the M31 
disk. First, we checked that the apparent magnitude of 
the BHB was consistent across all fields. While there 
was some increased scatter in F475W, it was small (<0.1 
mag), and the scatter was even smaller in F814W (<0.03 
mag). In neither case was a trend with galactocentric 
distance seen. Furthermore, we inspected the PHAT in- 
frared photometry of all of the regions where we detected 
the BHB, and found the red giant branch to have a single, 
narrow tip in all regions. We also checked the redden- 
ing value for the nearest globular cluster to the region, 
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which is very low (B201; E (B _ V) =0.04±0.02; ^MLStMi 
I2008T ) . Finally, wc inspected the Spitzer 24/x maps of the 
regions, finding no significant emission above the back- 
ground. Therefore, with all of these tracers, there is no 
evidence for the presence of a significant dust layer at 
these radii along the minor axis that would affect our 
measurements of the BHB. 

2.1.2. BHB Upper Limits 

Outside of 2.5 kpc on the minor axis and 5 kpc on the 
major axis, artificial star tests show that our 50% com- 
pleteness limit is a magnitude or more below the BHB 
feature, making it simple to isolate and measure the BHB 
component. In contrast, we were not able to isolate the 
feature inside of 2.5 kpc on the minor axis and 5 kpc on 
the major axis, because our photometric uncertainties 
and completeness affected our data too strongly at the 
magnitudes of interest. Attempts to solve this problem 
by applying completeness corrections from false star tests 
were not effective, due to the sharp cutoff in complete- 
ness in our region of interest. Essentially, when com- 
pleteness corrections become greater than a factor of 2, 
the inferred number of stars has associated uncertainties 
that our fitting routine could not overcome. Therefore, 
at these inner radii, we assume that the BHB is located 
in the same color-magnitude range as in all of the fields 
where it was measured directly (25.1<F814W<25.3 and 
0.1<F475W-F814W<0.5). We then counted the total 
number of stars in this range and corrected for the mean 
completeness. This calculation gives us the total number 
of stars which potentially could be BHB stars. However, 
because we cannot reliably subtract off the contribution 
from the MS or post-AGB stars, we quote our measure- 
ments in these fields as upper limits. 

Our interpretation of the inner radius measurements 
as upper limits is unaffected by magnitude biases from 
crowding. In very crowded fields, such as the inner M31 
bulge, our photometry shows biasing towards brighter 
magnitudes. Such biasing would put the BHB stars 
brightward of the magnitude range within which we 
are counting BHB stars. Because the luminosity func- 
tion increases toward fainter magnitudes, measuring the 
completeness-corrected numbers of stars at a brighter 
magnitude interval (to take such biasing into account) 
would only serve to lower our upper limits, making 
them less conservative. In one case, our completeness- 
corrected LF showed a small peak at F814W=25.05, 
slightly brightward of our assumed BHB magnitude. We 
verified that shifting our window 0.1 mag brightward 
(25.0<F814W<25.2) to include this small peak did not 
significantly impact our results. Such a change only 
slightly shifted our upper-limits, and the lowest upper 
limit changed by <3%. Thus our upper limits are con- 
servative and are not strongly sensitive to the precise 
magnitude cuts. 

We attempted to augment our radial baseline by ap- 
plying the same measurement techniques to other fields. 
We first attempted to detect the BHB in the PHAT sur- 
vey data along the major axis. However, we were unable 
to isolate a BHB feature because the expected number of 
BHB stars from the halo is small compared to the num- 
ber of main-sequence stars present in the disk component 
along the major axis. These fields had such high contam- 
ination levels that upper-limits measured using our blind 



counting technique were so large that they proved to be 
of no use. We therefore do not report any measurements 
taken near the major axis. 

2.2. Relevant Archival Data 

In order to greatly increase our radial baseline, we 
searched for relevant archival data that would contain re- 
solved photometry of the BHB in the M31 halo. First, we 
found the deep photometry of the halo released as high 
level sc ience products (PIDs: 9 453, 10265, 10816, PI: 
Brown. iBrown et al.l [20091 12008L hereafter B08). Then 
we found deep archival imaging away from the major 
axis or known streams that we could quickly process with 
the P HAT pipeline (PID: 10394, PI: Tanvir. lTanvir et al.l 
[20121: PID: 11362, PI: Rich) 

2.2.1. Public Deep Halo Photometry 

We downloaded the very deep photometry catalogs of 
B08, from ACS fields at 11, 21, and 35 kpc out along 
the minor axis, as well as one located at R=20 kpc in a 
known stream. These fields are well away from the disk, 
contain stars at the location of the BHB, and have no 
main-sequence stars.. 

For these data, which use a different filter set and dif- 
ferent photometric system, we found that the feature was 
isolated enough that a comparable measure of the num- 
ber of BHB stars (leaving out the instability strip and 
the vertical extension) was obtained from counting all 
stars with -1.0<F606Wsta/ag-F814Wsta/ag<-0.6 
and 26.0<F814W ST a/ag<26.7. The difference in 
F814W magnitude corresponds simply to the difference 
in F814W zeropoint between VEGAMAG (25.5) and 
STMAG (26.8). The color range was chosen by look- 
ing at the CMDs to find the color corresponding to the 
BHB feature; however, it is similar to the expected shift 
of —1 (a formal transformation of F475W-F814W=0.3 
results in F606W S tmag-FSUW S tmag= ~ 0.7). 

Although some of these fields are known to lie in 
streams, the fr action of stream sta rs is well-constrained 
by kinematics ([Gilbert et al.H2009b[ K For the 11 kpc field 
and 20 kpc stream field, only 56±16% and 25±10%ofthe 
stars b elong to the kinematically hot halo ([Gilbert et al.l 
l2009bl ). These fractions are determined by spectroscopy 
of RGB stars in fields that overlap the HST imaging. 
The kinematics of stars in these fields display distinct 
cold streams with low dispersion as well as a popula- 
tion of stars with high dispersion (hot) halo component; 
maximum-likelihood fits to the line-of-sight velocity dis- 
tributions in each field yield estimates of the fractions 
of stars in the hot and cold components. While both of 
these fields are in known streams, the outer field lies on 
a more overdense stream, resulting in a lower fraction of 
hot halo stars. Therefore, we scale our measurements in 
these fields by 0.56 and 0.25 respectively, to account for 
the overdensities, and we include the associated errors 
in our uncertainties. We note that accounting for these 
overdensities is necessary because of the small sample 
area of the ACS field of view at the distance of M31, as 
well as the bias of HST programs to point at streams. 
Larger areas in random halo locations, as one would ob- 
tain in more distant systems, would include both over- 
dense and underdense regions, averaging out to the true 
mean density of the halo. In addition, there is some 
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chance that the streams and the rest of the halo have 
different BH B properties in M31, as has been seen in the 
Galaxy (e.g. lBell et a,l.|[20Toh ; however, the similar Fbhb 
values seen through the streams and halo (see Figure QJ 
suggests that, if such differences are present in M31, they 
do not significantly affect our study. 

2.2.2. Other Relevant Archival Observations 

We searched for other archival ACS observations suit- 
able for our project. Although there have been many 
ACS observations in the M31 halo, most were either lo- 
cated on known overdense streams or were too close to 
the major axis to provide a clean measurement of the 
BHB. Unlike the B08 fields, other archival stream fields 
did not have spectroscopically determined stream frac- 
tions. Therefore we only included fields outside of known 
streams. 

We found 2 fields with favorable locations and depth. 
These were a parallel observation in program ff632, 
which was looking at background QSOs (we refer to this 
field as 11632JV131-HALO-SE), and an observation of a 
halo globular cluster (we refer to this field as 10394JV131- 
HALO-NW). For this latter field, a 20"x20" region cen- 
tered on the globular cluster we exclude from our analy- 
sis. 

These 2 fields were taken in F606W and F814W like the 
B08 fields; however, we reduced them in our own pipeline 
using a stricter quality cut appropriate for less crowded 
fields. In crowded fields, the crowding parameter is high 
for most of the stars because they all have close neigh- 
bors. In sparse fields, crowding is very sensitive to back- 
ground galaxies because shredded background galaxies 
have photometry similar to crowded stars. This con- 
taminant is insignificant in crowding-limited fields with 
hundreds of thousands of real stars, but is significant in 
sparse fields with tens of thousands of real stars. Much 
cleaner CMD features (such as the BHB) are gained by 
limiting the crowding cut in sparse fields. We produce 
CMDs similar in quality to those of B08 at the BHB by 
applying a crowding cut of O.f to these halo fields. Our 
artificial star tests showed that even with this more con- 
servative cut, we are complete at the location of the BHB 
in these sparse fields. 

3. RESULTS 

3.1. BHB Magnitude 

Our measurements of the BHB feature are detailed 
in Table [TJ These include, for each ACS field, the 
median (projected) galactocentric distance of the stars 
from the center of M31, the apparent magnitude of the 
BHB peak in our color range, the number of BHB stars, 
the number of RGB stars in our selection region, and 
the resulting BHB/RGB fraction. The measurements 
are plotted in Figures [3HU We find that the feature 
has an apparent magnitude of F814W=25.23±0.03, at a 
F475W— F81 4W color of 0.3. Assuming a distance mod- 
ulus of 24.47 (jStanek fc Garnavichlll998D and foreground 
extinction of Ay=0-21 ( Schlegel et al.lll998| ). this corre- 
sponds to Mj?gi4n/=0.63±0.05. The magnitude and color 
are consistent with the model BHB shown in Figure [TJ 
which has [Fe/H]= — 1.7; however, all model BHBs with 
— 2.3<[Fe/ H]<— 1.0 are similar to our observed BHB. 

3.2. BHB/RGB Fraction 



There is also information contained in the relative 
number of stars in the BHB feature compared to more 
well-populated features. For example, since BHB stars 
are typically produced by low metallicity populations, 
there may be a relation between the fraction of BHB 
stars present in an old population and the population's 
metallicity. To assess the relative strength of the BHB 
feature in our data, we calculated the ratio of the number 
of stars in the Gaussian component of the histogram fit 
(Nbhb) to the number of stars with 22.0<F814W<22.5 
and 1.5<F475W-F814W<3.5 {N RGB ). We took the 
RGB sample from 23.8<F814W S tmag <23.3 in the B08 
data. This magnitude slice provides a good proxy to the 
total stellar mass in the field, as it is dominated by the 
very well populated RGB. We tailored the magnitude 
range to yield a fraction (Fbhb = Nbhb/Nrgb) close 
to unity for ~12 Gyr old clusters with [Fc/H]~— 1. These 
values allow for easy comparisons between the halo light 
fraction and Fbhb (see § 13. 5p 

Once we had defined F bh b to measure the strength of 
the BHB in our data, we checked its sensitivity to pop- 
ulation metallicity by measuring it for a large sample of 
well-studied globular clusters. This comparison was per- 
formed to look for similarity between the M31 halo and 
Galactic globular clusters of similar age and metallicity. 
Since the Galactic halo, and most of the globular clusters, 
significantly differ in age and metallicity from the M31 
halo, only the higher metallicity and somewhat younger 
end of the globular clusters is directly comparable. How- 
ever, we i nclude the full ACS g l obular cluster trea sury 
program (jSaraiedini et all 120071; iDotter et etlj 120101) for 
context. 

Specifically, we measured the BHB fractions of the 
Galactic globular clusters (GC s) of the ACS glob- 
ular cluster treas ury program (jSarajedini et al.l 120071 : 
IDotter et al1l2010l ) All GC CMDs were corrected for their 
different dist a nces a nd extinctions using the values from 
IDotter et al.l (|2010f ) . Valu es for metallicity and ag e for 
NGC 6388 were tak en fromlWorlev fc Cottreil poTol ) and 
for NGC 6441 from lGratton et al l(|2006l ). Distances and 
redden i ng va lues for these two clusters were taken from 
iHarrisK) 19961 ). We then took our BHB from the CMD re- 
gion with -0.1<My-Mj<0.3 and -0.3<Mj<0.8, and 
our RGB from the CMD region with -2.0<Mj<-1.5. 
The results are shown in Figure [5j Applying the mean 
metallicities and ages from the full-field CMD analysis of 
B08, we find that the BHB fraction of the M31 halo is 
in good agreement with the FBHB-[Fe/H) correlation for 
Galactic GCs, despite the broad metallicity spread con- 
tained in the M31 fields^ Fbhb values for the M31 halo 
(~0.7) are similar to those measured for GCs of similar 
ages (11 Gyr) and metallicities (—1.0 0.5). Further- 
more, we note that Fbhb is the same in fields inside 
and outside of streams, suggesting that the streams have 
similar Fbhb to the kinematically hot halo. This find- 
ing is consistent with the full-field analyses of B08, which 
found the ages and metallicities of the steam field to be 
similar to other halo fields. 

Figure [5] confirms that at old ( ^ 10 Gyr) ages strong 
BHB populations are associated with low metallicities. 

9 Although the B08 metallicities are not from spectroscopy, they 
arc in agreement with spectroscopic metallicity measurements of 
stars in the vicinity. 
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Below [Fe/H]< — 1.5, we find Fbhb is universally high 
(1<Fbhb<6) with a median of 3.2. At higher mctal- 
licities, the strength of the BHB population drops dra- 
matically, falling by a factor of 10 as the metallicity rises 
to [Fe/H]<— 0.5. Above this metallicity, the BHB is un- 
likely to exist at all. Finally, we note that this compari- 
son includes the moderately metal-rich globular clusters 
NGC 6388, NGC 6441, and Lynga 7, which exhibit blue 
horizontal branches, showing that these clusters do not 
significantly affect the overall behavior of Fbhb with age 
and metallicity. 

We can use the data in Figure [5] to place limits on the 
metallicity distribution of the M31 stellar halo. If we as- 
sume that all of the BHB stars come from a metal-poor 
tail in the population and that the outer fields from the 
archive are pure M31 halo, we can infer from Figure [5] 
that no more than ^10% of the halo population is likely 
to belong to a metal-poor than with [Fc/H]< — 1.5. Such 
a tail would produce Fbhb values higher by a factor 
of 2 in these pure halo fields. This result is also con- 
sistent with the full CMD analyses summarized in B08, 
which have little stellar mass at [Fe/H]< — 1.5. We note 
that this result relies on the assumption that Fbhb be- 
haves in a similar way in M31 and the Galaxy. Figured] 
is consistent with such an assumption; how ever, there 
have been studies that suggest otherwise (e.g. lRich et al.l 
[2005ft . 

Finally, Figure [5] shows us that the youngest (and most 
metal-rich) halo field appears to have a somewhat high 
Fbhb value for its metallicity, making it consistent with 
the values from the other fields. Thus, our measurements 
suggest that the Fbhb value in the halo is relatively in- 
dependent of radius. Therefore, if the halo indeed has a 
gradient in its mean metallicity, as has b e en suggested 
in se veral works (e.g. IKalirai et al.l 120061 : IBrown et al.l 
120081) . then the gradient does not appear to strongly af- 
fect Fbhb in the radial range considered here. This re- 
sult is consistent with the fact that the BHB surface den- 
sity measurements obey a single-sloped power-law profile 
from 3-35 kpc in Figure UJ 

3.3. BHB Surface Density Profile 

We now analyze the surface density profile of BHB 
stars. If the stellar halo can be defined as the population 
that is responsible for the field BHB stars, then this pro- 
file should be consistent with previous measurements of 
the outer halo. We show our measurement of the surface 
density profile of BHB stars in M31 in the left panel of 
Figure SJ The surface density falls off with galactocentric 
distance following roughly the surface brightness profile 
of the halo measured by G05 and Cll outside of 11 kpc. 
However, the BHB profile is steeper between 3 kpc and 
11 kpc than inward extrapolations of such current halo 
models. Interestingly, our data suggest a single power- 
law slope all the way from 3 kpc to 35 kpc. The slope is 
steeper than (but within the uncertainties of) G05, and 
it is similar in slope to the outer halo portion of the fits 
of Cll. Thus, our data are providing a new constraint on 
the shape of the inner halo, suggesting the single power- 
law slope extends inward to 3 kpc before flattening. 

We fit to the BHB surface density profile using the 
halo parametrization of Cll. When fitting the data, we 
only applied our lowest measured upper limit by setting 
the value and uncertainty both equal to half of the up- 



per limit. The power- law function from Cll provides an 
excellent fit to the BHB data (x 2 =0.97). 



E(r) 



(l+r/a h Y 



a/2 



Where E(r) is the BHB surface density at radius r, a is 
the projected 2D spatial density distribution power-law 
index, and ah is a core radius (in units of kpc) inside 
of which the profile flattens. E* and -R* arc normal- 
ization constants. Not surprisingly, the best fit has the 
surface density right at our lowest upper limit. Our data 
constrain the model to a=2.6to 2 an d flft,=2-7^J'g kpc. 
While this functional form features a core and a central 
slope that goes to zero as r — > 0, we do not claim the 
actual detection of a core in the inner halo, only that 
the slope decreases to <1.2 interior to 3 kpc ( a steeper 
profile would exceed our upper limit at 1.6 kpc). We 
only use this functional form to aid in comparison with 
previous measurements. 

How do our halo parameters compare with other mea- 
surements? We first compare with other estimates of 
power law slopes. Cll's surface brightness decompo- 
sition study gives a power law index of a=2.52±0.08, 
in excellent agreement with our power law. G05 finds 
a = r~j — 2 with no quantitative uncertainties, and thus 
roughly consistent with our result. A more formal result 
from a spectroscopically confirmed sample of M31 halo 
stars yields a = 2.2±0.2 (K. Gilbert et al., ApJ, sub- 
mitted), also consistent with our measurement. We have 
also estimated a flattening radius for the BHB profile (a 
robust measurement), which as long as the BHB/RGB 
ratio for the halo is constant (an assumption) is consis- 
tent with a core radius for the stellar halo as a whole. 
Our estimate does not agree with others in the litera- 
ture (Cll have a large core, attributing instead much of 
that light to the disk component; and G05 lack a core 
at all in their parametrized profile), but the constraints 
from other works were necessarily weak due to the small 
fraction of halo stars at these small radii. 

3.4. The Mass of M31 's Stellar Halo 

We can integrate our 2-D projected surface density pro- 
file to obtain an estimate of the total halo stellar mass. 
First, we make the strong assumption that the BHB sur- 
face density is directly proportional to the halo stellar 
mass surface density, effectively making E(r) the profile 
of the stellar halo surface density. Taking /xy=29 mag 
arcsec -2 at 21 kpc (from integrating the 21 kpc field 
catalog of B08), My =4.82 and assuming a stellar mass 
to V-band light ratio of 1.5 dBcllazzini et al.l I2012T) . we 
calculate E*=1.4xl0 5 M Q kpc" 2 and i?*=21 kpc. 

Then, we assume azimuthal symmetry, making the to- 
tal mass equal to the following integral: 



/^2tt ,-260 
E(r)dA = / / rE 
Jo Jo 



(r) droie 



This calculation yields a total halo stellar mass out 
to the virial radius (260 kpc iSeigar et al.l |2008|) of 
2.ll^xl0 9 Mq, with a half- mass radius of 7 kpc. 
Changing the outer boundary of the integration has 
a small effect on the resulting value. Integrating to 
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40 kpc (near the outer boundary of our photometry data) 
yields 1.8 xlO 9 M Q , while integrating to infinity yields 
2.3 xlO 9 M©. We also note that many other stellar ha- 
los/envelopes appear to be oblate within 10 — 20 kpc 
with a xis ratios c/a betw een 0.4 and 0.7 (e.g., the Milky 



Uuric et al.l 



Barker et al 



Barker et al 



20081: NGC 253: IBailin et al.l [20TTI: 
2009, M. Vlajic et al. in prep.; NGC 



2012), it is possible that our assump- 



Way: 
M81: 
2403: 

tion of azimuthal symmetry will turn out to be incorrect. 
In that case, we expect that the stellar mass estimate 
would scale very approximately as a/c (the inverse of 
the projected axis ratio of the stellar halo). This sim- 
ple scaling applies to the idealized case, where the M31 
minor axis is aligned with that of the minor axis of the 
spheroid and there are no projection effects. 

3.4.1. Comparing with Other Stellar Halo Measurements 

Comparing with previous estimates of M31's halo 
mass, our measurement is slightly on the high side, but 
broadly consistent with expectations. We note that our 
estimate is less than an inward extrapolation of the 
power-law (index = 2.2±0.2) of Gilbert et al. (ApJ, sub- 
mitted), or other s imilar results (e.g. index=2 . 1 7±0 . 1 5 , 
iTanaka et aT1l20Tot index=1 91±0. 11, llbata et al.ll2007l : 
index~-2.3. llrwin et al.ll2005h which results in a halo stel- 
lar mass (assuming no core and excluding the cusp inside 
of 1 parsec) of ~1.3xl0 10 Mq out to the virial radius 
and a half- mass radius of <1 kpc. If we apply the ah 
value measured from the BHB profile to the Gilbert et 
al. (ApJ, submitted) profile, the resulting halo stellar 
mass will be the consistent with our value (1.8xl0 9 Mq, 
but with a half- mass radius of ~16 kpc ) . As suming 
M/L~2, the measuremen t of llbata et al.l (|2007l) yields 
-2 xlO 9 Mq. Ilrwin et al.l (|2005| ) estimate ~2.5% of the 
total light from M31 is from the halo, and Cll estimates 
4%. Assuming a total M31 luminosity of 2xl0 10 L BQ 
and M/Lb~2, these correspond to stellar halo masses of 
1.0-1.6 xlO 9 M Q , respectively. 

We can also compare our M31 stellar halo mass with 
those of other nearby large disk galaxies. Our estimate 
of the stellar mass of the halo of M31 is less than cur- 
rent estimates for the mass of the stellar halo of the 
nearb y disk galaxy NGC 253 (4xl0 9 M Q IBailin et al.l 



120 111) , but larger than current estimates for the mass 
of the stellar halo of the Mil k y Wa y out to 40 kpc 
(3.7±1.2xl0 8 M ? . iBell et all 12001 2-10xl0 8 M r 
iSiegel et al.l[200l Uuric et al.ll2008t iDeason et al.ll20Tll 

3.4.2. Comparing with Simulations 

M31 and the Milky Way live in massive dark matte r 
halos of similar mass (1.4xl0 12 Mg lWatkins et al.ll20l"of ). 
The stellar mass of the Milky Way disk is estimated to 
be 5.5xl0 10 M Q (jFlvnn et alJl2006l ). and that of M31 is 
estimated to be roughly the same when deter mined from 
the a bsolute visual magnitude of My~- 21.0 (jFont et al.l 
1 2 1 If) . Thus, simulations of Milky Way-like galaxies can 
be compared to M31 observations as well as Milky Way 
observations. 

Comparing to cosmologically-motivated simulations in 
which the stellar halo is built up through the accretion 
of dwarf galaxies alone, we find a reasonable degree of 
agreement between the properties we infer for M31's 
stellar halo and BJ05 and C10. BJ05 find halo stel- 



lar masses ~ 2 x 10 9 M Q , with 3-D power law slopes 
outside lOkpc of between —2.5 and —3.5. C10 find 
a somewhat wider range in possible halo masses, be- 
tween 10 8 M Q and 2.5 x 1O 9 M , and 3-D power law 
slopes ranging between —2 and —4.5. Hydrodynami- 
cal models of stellar halo formation that include an in 
situ componen t find rather higher stellar masses, 0.5- 
3.4xl0 10 M Q (jZolotov et al.l [2009IF I. Possibly relevant 
given the strong evidence of a very active M31 merger his 



tory fe.g. iGilbert et aj . 2009a; McC onnachie et a 



is that both IZolotov et al.l (j2009h and iFont et al 



2009) 



20111 ) 



find that the galaxies with the most active merger histo- 
ries have the lowest in situ fractions in their simulations. 
Yet, simulated stellar halos have large halo-to- halo scat- 
ter which is merger history-dependent; given this scatter, 
we conclude at this stage that the properties of M31's 
stellar halo appear to be consistent with the range of 
stellar halos expected in a cosmological context. 

3.5. BHB/RGB Ratio Profile 

We can use the data in the right panel of Figure |4] 
to assess the fraction of halo stars at small radii. In 
the right panel of Figure 0] we plot Fbhb as a function 
of galactocentric distance, along with the G05 and Cll 
models of the fraction of total light coming from the M31 
halo. The models suggest that the relative strength of 
the halo clearly increases with galactocentric distance out 
to 11 kpc, then remains relatively constant at the pure 
halo value. This radial distribution suggests that the 
population becomes dominated by halo members outside 
of ~11 kpc. 

The behavior of Fbhb in Figure |4] is consistent with 
what we expect for a halo population based on the G05 
and Cll models. At large radii, the measured values of 
Fbhb are consistent with being flat (although the un- 
certainties are large, due to the small number of stars). 
At smaller radii, however, the strength of the BHB fea- 
ture falls dramatically compared to the RGB. Empiri- 
cally, this drop indicates the increasing contribution of 
old stellar populations that do not host BHB stars. If 
Fbhb for the M31 halo is roughly constant with galac- 
tocentric distance, then the low values inside of 11 kpc 
are due to the increased RGB contribution from the bulge 
and disk components. If we assume that the BHB stars 
are all halo members and that Fbhb for the halo is 0.7, 
then the expected behavior of Fbhb is to follow 0.7 H/T 
where H is the halo stellar surface density and T is total 
stellar surface density. 

Perhaps not surprisingly, both the radial distribution 
of Fbhb and the number density profile suggest that the 
BHB stars inside of 4.5 kpc do not fit inward extrapola- 
tions of current models of the halo profile. Inside of 5 kpc, 
Fbhb is significantly larger than expected. This excess 
in BHB stars suggests a number of possibilities, which we 
have explored, or address in the subsections that follow. 
1) The halo has a different shape from inward extrapo- 
lations of current decompositions (the interpretation we 
explored in sections 3.3 and 3.4). 2) Fbhb of the halo 
is not constant. 3) The bulge is contaminating the BHB 
sample. 4) The disk is contaminating the BHB sample. 
We now discuss these latter 3 possibilities and argue that 

10 IFont et al.l Q201H ) do not quote stellar halo masses as all the 
properties quoted in their work are for bulgc+halo. 
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they arc unlikely, leaving us to conclude that we have in- 
deed identified a means by which to constrain the shape 
of the inner stellar halo. 

3.6. Does the Halo Produce BHB Stars at all Radii? 

An assumption involved in our analysis is that the M31 
stellar halo produces field BHB stars in the same fraction 
at all radii. While our results from 10-35 kpc suggest 
that this assumption is reasonable, there is no clear way 
to test the assumption at small radii. Indeed, our upper- 
limits inside of 2 kpc allow the possibility that the there 
is no population that produces a field BHB at these small 
radii. We stress that our profile only describes the popu- 
lation that produces a field BHB, which we suggest may 
be a reasonable way to define the M31 stellar halo. Thus, 
in this analysis, all kinematically hot stars beyond those 
required to produce the BHB profile shape (i.e. many of 
the stars inside of 3 kpc) are associated with the M31 
bulge. 

Even if the halo produces Field BHB stars at all radii, 
Fbhb of the halo component may not be constant. The 
outer M31 halo (outs ide of ~10 kpc) is known to have a 
metallicity gradient (|Kalirai et al.l [2006h . However, this 
gradient is negative (metallicity increases with decreas- 
ing radius) . While it is known that the total HB fraction 
(inc luding the RHB) i ncreases with increasing metallic- 
ity (jSalaris et al.ll2004l) . based on the behavior of Fbhb 
in Figure |U Fbhb decreases with increasing metallicity. 
Thus, we would expect the stellar halo only Fbhb to de- 
crease slightly toward the galaxy center. Such behavior 
would make Fbhb even lower than the dashed or dotted 
curves in Figure 4 near the center, which is the opposite 
of the observed behavior. Therefore, a changing Fbhb 
of a halo of the type currently in the literature appears 
an unlikely explanation for the observed radial profile. 
Wc note that if the halo Fbhb does decrease towards 
smaller radii, the stellar mass and power law slopes of 
the halo inferred in sections 3.3 and 3.4 would increase. 
We now move on to discuss the possibility of bulge or 
disk contamination. 

3.7. Are the BHB Stars From the Stellar Halo? 

An important assumption involved in our use of the 
BHB as a probe of the M31 stellar halo is that neither the 
bulge nor the disk components contribute significantly to 
the field BHB. Here we argue that our data suggest this 
assumption is reasonable. As detailed below, if we as- 
sume that a significant fraction of BHB stars belongs 
to either the bulge or the disk, we can fit the profile, 
but only if these components contain old, metal-poor 
populations that have not been observed any other way. 
Furthermore, the density of BHB stars from 3-35 kpc 
appears to follow a single power-law distribution. Thus, 
the BHB stars appear to belong to the stellar halo, which 
extends to small radii with a slightly different structure 
than inward extrapolations of currently-available decom- 
position models. Inside of ~3 kpc, simple inward ex- 
trapolations greatly overpredict the observed numbers of 
BHB stars, suggesting a break in the profile. 

3.7.1. Possible Contributions of Bulge BHB stars 

One possible source of excess BHB stars at small galac- 
tocentric distances is the bulge. However, given the 



degree to which the bulge dominates at small radii, 
the low numbers of BHB stars in M31 seen inside of 
2 kpc suggests that the bulge population is not respon- 
sible for the excess. Furthermore, the high- metallicity of 
the b ulge (0.0 < [Fe/H] < 0.4 out to 1 kpc, iSaglia et al.l 
l2010h also argues against it as the source of the BHB 
stars seen at 2.7-4.5 kpc. On the other hand, the 
bulge is known to contain a b road range of metallicities 
(ISaraiedini fc Jablonkal[2005h . If we assume that some 
fraction of the bulge population is metal-poor enough to 
produce a significant BHB ([Fe/H]<— 0.5), then we can 
estimate what fraction of metal-poor stars would best 
reproduce the BHB profile. If the resulting estimate is 
reasonable, then our assumption that the BHB traces the 
halo becomes questionable. 

The potential contribution of the bulge to the BHB is 
different depending on the M31 bulge profile. Therefore 
we tested both the bulge profiles of Cll and G05, which 
bracket the possibilities of halo size. The Cll bulge is 
small, and the G05 bulge is large. 

First, assuming the decomposition models of Cll, we 
can calculate the bulge properties necessary to reproduce 
the observed BHB density profile. The preferred Cll fit 
has a small bulge component that decreases very steeply 
to zero by a minor axis distance of 4.5 kpc. Given this 
decomposition fit, it is impossible for the bulge to be 
responsible for the BHB stars because there would not 
be enough stellar mass associated with the bulge at 3-5 
kpc to produce the observed number of BHB stars, even 
if the bulge metallicity were as low as the halo metallicity. 

If instead we assume the decomposition model of G05 
(large bulge), we can produce a BHB profile consistent 
with the observations if the bulge has a steep increase 
in the fraction of metal-poor stars outside of 1 kpc. If 
we assume that some fraction of the bulge population 
is metal-poor enough to produce a BHB ( [Fe/H] <— 0.5), 
we can estimate what fraction of metal-poor stars would 
best reproduce the BHB profile. Our very conservative 
upper-limits at these radii (1.6 kpc) suggest that less that 
10% of the inner bulge population is metal-poor enough 
to produce a BHB. We can bring the profiles in complete 
agreement by imposing a fraction of metal-poor stars in 
the M31 bulge that increases very steeply, from 0% to 
50% from the galaxy center to 2.5 kpc. Another way 
to bring the profiles into agreement would be to decrease 
the metallicity of the metal-poor component of the bulge. 
This explanation is possible, but it is not likely, as it 
would be odd if the metal-poor stars of the bulge were 
of lower metallicity than those of the halo. Since no 
such steep metallicity gradient has been observed at these 
radii in the bulge, it appears unlikely that the BHB stars 
at 2.7-4.5 kpc are dominated by the bulge component. 

3.7.2. Possible Contributions of Disk BHB stars 

Another possible source of excess BHB stars is the 
outer disk component. However, there is a very powerful 
empirical argument against disk stars being responsible 
for the prominent BHB population at 2.7—4.5 kpc along 
the minor axis: the BHB is not detected in the outer 
fields along the major axis of the PHAT survey data, as 
would be expected if the disk stellar populations are well- 
mixed (as an old population like the BHB is expected to 
be). The deprojectcd disk radius of 2.7-4.5 kpc is 1-1.5 
degrees. The PHAT data only extend to ~1 degree along 
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the major axis, but the B08 disk field (~2 degrees out 
on the major axis) shows very little evidence for the old 
(>10 Gyr), metal-poor component (| Brown "et~aT1[2006h 
needed to produce a BHB excess. Therefore it appears 
that the disk is unlikely to be a significant source of BHB 
stars. It is worth noting that this lack of BHB stars is 
not necessarily because the metallicity of the outer disk 
isn't low enough, but rather may indicate that the vast 
majority of outer disk stars are too young to produce a 
prominent BHB. Furthermore, kinematic measurements 
suggest a low disk c ontribution at these minor axis radii 
(jGilbert et al.l [20071 < 10% at 9 kpc). Thus, both the 
observed populations of the disk itself and the spectro- 
scopic disk fraction argue against a disk origin of the 
BHB stars at these radii. 

4. CONCLUSIONS 

We have measured the number density of BHB stars in 
M31 at galactocentric distances ranging from 1.6 kpc to 
35 kpc using photometry from the PHAT survey along 
with archival halo ACS fields. Our measurements show 
that the properties of the BHB of M31 are consistent 
across two degrees of galactocentric distance. Galactic 
globular clusters with ages of ~10 Gyr and metallicities 
of -1.0<[Fe/H]<-0.5 (similar to those of the M31 halo 
values measured from deep CMD analysis) have BHB 
properties that match the M31 halo as well. Taken to- 
gether, these measurements suggest the BHB is a useful 
tool for tracing galaxy halos. 

In additional tests, we showed that the BHB surface 
density profile follows that of the known halo outside of 
10 kpc. However, there is an excess of BHB stars at 2.7- 
4.5 kpc over inward extrapolations of current M31 halo 
profiles. To match the radial profile of the BHB a combi- 
nation of halo and bulge and/or disk components can be 
applied; however, the contributions of these components 
that is required is inconsistent with the expected prop- 
erties of the M31 bulge and disk populations at these 
radii. While it is possible that the disk and/or bulge 
populations could be contaminating our BHB sample, 
there is no compelling evidence that the M31 bulge or 



disk harbors a significant BHB population at any radius. 
Furthermore, these components would need to be con- 
tributing BHB stars in just the right proportions for the 
BHB density to follow a power-law density distribution. 
Such a conspiracy is not impossible, but seems unlikely. 
Thus, our use of the BHB to trace the inner halo is fur- 
ther justified. 

The BHB data are well matched by a power-law with 
an index of 2.6+qJ outside of 3 kpc decreasing to <1.2 
inside of 3 kpc. This profile describes the population that 
produces a field BHB, which we suggest is a reasonable 
way to define the M31 stellar halo. In this picture, all 
kinematically hot stars that do not follow this profile 
could be associated with the M31 bulge. Our profile 
slope is consistent with the range in slopes characteristic 
of simulations of stellar halo formation in a cosmological 
context. 

Normalizing our best-fitting profile function to a halo 
stellar mass density of 1.4xl0 5 M Q kpc -2 at 21 kpc 
yields a total stellar halo mass of 2.1^4 xlO 9 M for 
M31. This mass compares well with, but is on the high 
side of, other estimates of the M31 halo stellar mass, 
and it is significantly higher than current estimates of 
the stellar halo mass of the Milky Way. 

About a decade ago, the fir st large stream in M31 was 
discovered (jlbata et alJl20Qll ). Since then, the M31 halo 
has been mapped to large radii with Keck (K. Gilbert et 
al., A pJ, submitted) and the CFHT (|McConnachie et al.l 
l2009f ) . and all measurements show significant structure 
on small scales indicative of a very rich merger history. 
The BHB stars of the halo as measured with HST ap- 
pear to paint the same picture. The shape and mass of 
the stellar halo as measured with BHB stars are broadly 
consistent with those of simulated halos that contain a 
significant fraction of accreted stars. 
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TABLE 1 

Properties of the M31 Blue Horizontal Branch 
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a The median galactocentric distance of the stars in the the region kpc. 

b The F814W magnitude of the BHB at 0.1<F475W-F814W<0.5. 

c The number of BHB stars from 0.1<F475W-F814W<0.5 or equivalent. 

d The number of RGB stars from 1.5<F475W-F814W<3.5 and 22.0<F814W<22.5 or equivalent. 
e Fraction of BHB/RGB stars. 
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12 3 

F475W - F81 4W (vegamag) 

Fig. 1. — The color-magnitude diagram of the first PHAT field where we noticed the BHB (R m inor=3.8 kpc). Redder colors denote 
higher densities of data points. Ovcrplotted in white is the theoretical zero age horizontal branch for [Fe/H]= -1.7 from the Padova stellar 
evolution models. The overdensity in the CMD associated with the BHB is clearly evident as the density peak coincident with the model 
at F475W-F814W~0.3, F814W~25.3. The metallicity has little effect in this sequence at the color of the BHB, but changes the shape of 
the red end. Ovcrplotted in dashed gray lines are the boundaries over which we fit histograms in order to measure the properties of the 
BHB stars. The dashed gray box outlines the area in which we counted RGB stars for BHB/RGB ratios. 
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Fig. 2. — The northern half of the GALEX NUV image of M31. Thin line boxes show the regions of the PHAT survey completed when 
the BHB was seen. Thick black box shows the region where the BHB was noticed in the photometry: the area farthest out on the minor 
axis, with the lowest disk contribution. 
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Fig. 3. — Left: Luminosity functions for stars with 0.1<F475W-F814W<0.5, in PHAT fields with clear BHB features. Overplotted are 
the best-fitting line+Gaussian functions. Each field is represented by a different color. The numbers in the upper left indicate the number 
of stars (and uncertainty) in the Gaussian component. Right: The resulting best-fit apparent F8f 4W magnitude for the center of the BHB 
as a function of galactocentric distance. Cyan points show the combined measurements of multiple fields at similar radii. There is no 
trend in BHB peak magnitude with radius or surface brightness in M3f , thus the Gaussian fits are not significantly affected by crowding or 
completeness, and that the detected stars in each location are behind equal amounts of extinction. The BHB in M3f has F814W=25.23. 




Fig. 4. — Left: The number density of BHB stars as a function of M3f galactocentric distance. Diamonds show combinations of the 
multiple measurements at similar radii in the PHAT data. Squares show the B08 and archival photometry. Upper-limits are shown with 
downward arrows. Dotted line: M31 halo profile from G05 normalized to match our data at 20 kpc. Dashed line: M31 halo profile from 
Cf 1 normalized to match our data at 20 kpc. Solid line: a fit of the BHB data using the power-law halo parametrization of Cf f (x^=0.97). 
Right: Fbhb as a function of galactocentric distance over the same radial range. Lines are from the same models as Left but represent 
0.7x H/T where H is the Halo light and T is the total light in the two models. The low BHB/RGB ratios near the center are due to the 
bulge and disk components contributing substantially to the RGB but not to the BH B, a nd do not reflect metallicity changes in the halo 
component. The factor of 0.7 scales the fraction to the p ure halo ¥~gjjg value (see § 13. 51 . The curve corresponding to the BHB fit (solid 
curve) assumes the total from Guhathakurta ct al. (2005) to compute fractions and is only computed out to the "pure halo" value (0.7) at 
21 kpc. 
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Fig. 5. — BHB fractions measured for Galactic GCs from the sample of lDotter et al.l H2010T ) using the same CMD regions as for the M31 
halo fields of B08. Points are color-coded by age, given by the color bar in Gyr. GCs are shown with asterisks. M31 points are shown 
with open squares using the ages and metallicities measured in B08. Dashed vertical lines mark the metallicity of the M31 thick disk 
llCollins ct al. 201l|). Dashed horizontal lines lines mark the range observed at 3—5 kpc from the galaxy center. We assume that the low 
values at these inner radii are in large part due to the disk and bulge components adding to the RGB, but not to the BHB. 



